Critical Temperature for Bose-Einstein condensation in quartic potentials 
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The quartic confining potential has emerged as a key ingredient to obtain fast rotating vortices 
in BEC as well as observation of quantum phase transitions in optical lattices. We calculate the 
critical temperature Tc of bosons at which normal to BEC transition occurs for the quartic confining 
potential. Further more, we evaluate the effect of finite particle number on Tc and find that /^Tc/Tc 
is larger in quartic potential as compared to quadratic potential for number of particles < 10^. 
Interestingly, the situation is reversed if the number of particles is > 10^. 
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I. INTRODUCTION 



Particles in anharmonic potentials is a well studied ex- 
ample of chaotic system. A simple but good model of a 
system which exhibits parametric transitions from regu- 
lar to chaos and vice-versa is the quartic oscillator. The 
parameters which control the transitions are the coupling 
strengths of the cross terms. The quartic oscillator be- 
sides being a case study of chaotic system, plays a very 
important role in the guise of Higg's field, through which 
all the known fundamental particles acquire finite masses. 
The same potential also appears in optical lattices, where 
quantum phase transition from superfluid to Mott insu- 
lator has been observed Optical lattices are regular 
intensity patterns of light created with counter propagat- 
ing laser beams. The laser beams have a Gaussian profile 
and create an overlapping anharmonic potential across 
the optical lattice. Recent theoretical studies show that 
realization of the quantum phase transition can be more 
robust with a quartic potential ^ . Besides the quantum 
phase transitions, introducing a quartic confining poten- 
tial stabilizes fast rotating vortices in Bose-Einstein con- 
densates d, 13, d, H, 0j 9 ■ recent work [o^ , the con- 
densation temperature and thermodynamic properties of 
a rotating ideal Bose gas in an anharmonic trap has been 
studied. However, in the experimental realization of the 
quadratic plus quartic confining potential, the observa- 
tion of the fast rotating vortices eludes an unambiguous 
detection 

In this paper we calculate the critical temperature Tc 
at which bosons confined by a quartic potential condense. 
BEC in such a potential was studied for low-dimensional 
systems [ll|. This calculation requires the density of 
states. In our calculations we use the semiclassical ex- 
pression of the density of states, which is valid at higher 
energies. It is however not appropriate to study low lying 
states. In particular, the energy of the ground state is 
essential to estimate Tc for finite number of particles. To 



estimate the correction to Tc in the finite particle case, 
we calculate the ground state energy analytically [13] ■ 
The calculation is based on a method which optimizes 
the matrix elements of the quartic potential Hamiltonian 
in the harmonic oscillator basis. 

Our calculations show that Tc in quartic confining po- 
tential is higher than quadratic potential. This is per- 
haps to be expected, since Tc varies as iV^/" when the 
density of states is proportional to e"~^. In case of 3D 
harmonic oscillator potential a = 3 whereas it is 9/4 in 
the case of 3D quartic potential. Hence, in the quartic 
potential case Tc varies as A^^/^ compared to N'^^^ in the 
case of harmonic potential. From our calculations, it is 
evident that the cross terms increase Tc in 2D as well 
as 3D quartic potential. We find that Tc rises by factor 
of 1.2 and 1.1 in the 3D and 2D potentials respectively. 
However, the experimental realizations of the trapping 
potentials, which are created from laser beams, are more 
appropriately described without the cross terms. 



II. Tc FOR 3D QUARTIC POTENTIAL TRAP 

The general form of the quartic oscillator potential is 
A(r.r)^. It is homogeneous and has cross terms in the 
Cartesian coordinate representation which couple mo- 
tions along different axes. Potentials of this form occur in 
optical lattices, where counter propagating lasers create 
undulating patterns of standing radiation field. In one 
dimension, a pair of counter propagating Gaussian laser 
beams along z-axis of intensity profile /q exp(— 2r^/w) 
creates an array of periodic intensity minima and max- 
ima. These are located along the z-axis. Depending on 
the detuning of the laser, the atoms are attracted to the 
intensity minima or maxima. Usually the wavelength of 
the laser A is much smaller than the beam width w. To a 
very good approximation, the potential across a surface 
normal to the laser beam 
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Tuning the parameters of the laser beams, it possible to 
retain only the quartic term. For simplicity, neglecting 
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the cross terms, in three dimension 

V{x,y,z) = X{x^ + y^ + z^). 



(2) 



The eigen energies of the corresponding Hamiltonian is 
the sum of eigen values corresponding to each dimension. 
For the one dimensional quartic oscillator, eigen values 
can be calculated by minimizing the expectation of the 
Hamiltonian in the basis states of harmonic oscillator of 
appropriately chosen frequency p^. Generalizing the re- 
sult to three dimensional case, the eigen energy 



e(ni,n2,n3) = 1.389 VK + -)3 
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bmce 

1.389(A;i'*/m2)i/3 has the dimensions of energy we 
can represent this factor by huj, then 



e(ni, 712,713) 
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We now determine the number of states G{e) with energy 
less than a given value e. For energies large compared to 
Sw, we may treat t^^'s as continuous variables and neglect 
the ground state energy. To calculate G(e), we introduce 
a coordinate system in terms of the three variables — 
n^^^huj. In this coordinate system e — ei + e2+e3 defines 
a surface of constant energy e. Then G(e) is proportional 
to the volume in the first octant bounded by the surface 
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the integral we use the 



To evaluate 
/ X {u — 

B{ti,u) = r(/i)rH/r(M 

sity of states 



C-C1-C2 
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relation 
where 



[13]. Then the den- 



5(e) 
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de 



0.6852 



,5/4 
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This expression for density of states is used to calculate 
Tc- For bosons, the total number of particles occupying 
the excited states is 



iVn- 



5(e) 



;(e~M)/fcT _ I 



de. 
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At critical temperature /i ^ in the case of bosons and 
Nexc is equal to the total number of bosons N. Eval- 
uating the integral gives Tc in terms of the number of 
bosons 



kTr 



[0.6852r(9/4)C(9/4)f/'^' 
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The corresponding expression for three dimensional 
isotropic harmonic oscillator potential is 



(9) 




FIG. 1: Ratio of the critical temperature (Tc) to angular fre- 
quency (tJo) as a function of the number of bosons on the 
log scale. Here we have taken 1.389(A/l"'/m^)^''^ = hcuo {^^o 
is angular frequency of isotropic harmonic oscillator). When 
logAf — 7 , the Tc for the 3D quartic potential is approxi- 
mately six times higher than the 3D isotropic harmonic po- 
tential. 



From eq. ^ and eq. ^ , the ratio of the critical temper- 
atures in the two potentials is 



(-^c)quartic 



^4/9^(3)1/3 



armonic 



iVi/3 [0.6852r(9/4)C(9/4)f/^ 
= 1.0067V^/^. (10) 



Where we have taken huj — hwo to obtain the ratio. The 
ratio is proportional to N^^^, which means that for 10^ 
atoms, Tc in the case of 3D quartic potential trap is ap- 
proximately six times higher than that of the 3D isotropic 
harmonic trap. 



III. EFFECT OF FINITE PARTICLE NUMBER 

The expression of Tc in the previous section is with 
the approximation that the ground state energy is zero, 
which is a valid approximation when the system has large 
number of bosons. For finite number of bosons, the zero 
point energy causes a change in the value of Tc- For the 
3D quartic potential, the ground state energy is [l^ 



2.41 



1/3 



(11) 



This should be equal to the change in the chemical po- 
tential at the critical temperature, that is A/i = Cmin- As 
the total number of bosons is fixed [141 



dN = 



dN 
9T 



dT 



d^i J J, 



(12) 



3 



This implies 



df 
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Using the expression of N at temperatures shghtly above 



e(^-^)A^ - 1 



de. 
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This relation is obtained by substituting the general ex- 
pression for the density of states i.e. g{e) — Cqe""^ in 
eq. ([7]). Here Cq is a constant whose value depends on 
the form of the trapping potential. Then from eq. ([Te 
and eq. p4|) we get 



9/1 
df 



N 



C(«) , 

-a— -k. 

C(a-l) 



(15) 



In this expression a should be greater than 2, otherwise 
the relation is not valid since C(l) diverges. Using this 
expression, the change in the critical temperature due to 
the finite particle number is 



AT, 



C(«-i) 

aC(a)fc 



A/i. 



(16) 



In the case of 3D quartic potential a is equal to 9/4. 
Then 



AT, 



-4C(5/4) ^ 
9C(9/4)fc ""''^ 

-1.071C(5/4) /Aft' 



C(9/4)fc V m 



4\ 1/3 



(17) 



A relative measure of the effect of zero point energy on 
the critical temperature is the fractional change of the 
critical temperature. It is the ratio between ATc and Tc, 
for the present case 



AT, 



0.9054C(5/4) (Xh^V'"^ (r(9/4)C(9/4))*/^ 



C(9/4) 



1 /3 

Noting that {Xfi^/m'^^ is equivalent to fiu)/1.389 we 
get 

ATc -0.6891C(1.25)iV-4/9 



(18) 



C(2.25)5/9 
-2.56iV^'*/^. 



(19) 



(20) 



For the 3D isotropic harmonic potential, the fractional 
change of the critical temperature is 

AT, 



T, 



= ~0.73iV 



-1/3 



(21) 



If we compare eq. ([^U]) and eq. pT|) . we find that the per- 
centage decrease in T, is larger in the case of 3D quartic 
potential trap for number of particles ^ 80, 000. But the 
scenario is reversed for number of particles > 80, 000. 
This is also evident from fig. ^ where cross over point 
corresponds to the number of particles w 80, 000. 




FIG. 2; Fractional change in critical temperature (ATc/Tc) 
as a function the logarithm of the number of bosons (log A*" 
). The green and red colored plots correspond to the quartic 
and quadratic potentials respectively. Below log(A'^) = 5, that 
is A*' ~ 80, 000, the 3D quartic potential has larger fractional 
change. Above this point the 3D isotropic harmonic potential 
has marginally higher fractional change. 



IV. EFFECT OF CROSS TERMS 
A. 3D case 

Consider the general form of the quartic potential, as 
mentioned earlier 



V{r) = \{r.rf 



(22) 



In optical traps, it is possible to create confining poten- 
tials which are approximately close to this form but a 
truly spherically symmetric one is not realizable. The 
difhculty is in producing the cross terms of the potential, 
for example, terms like x^y^ in Cartesian coordinate rep- 
resentation. The absence and presence of the cross terms 
in quadratic and quartic potentials respectively introduce 
a key difference between the dynamics in the two poten- 
tials. In absence of the cross terms like xy, in quadratic 
potential, a perturbation to the dynamics of a particle 
along an axis remains confined along that axis. In con- 
trast, it propagates to other axes in the case of quartic 
potential. For condensates in traps, an important pa- 
rameter which reflects the effects of these terms is the 
critical temperature. 

Semiclassically, total number of states available to the 
system can be obtained by dividing the total phase space 
volume by /i^, the volume of a single state 



G(e) = 1 / dx / dp. 



(23) 



= ^ / ^^d^ / P^dp. (24) 
h'^ Jo Jo 

In the above equation r* and p* are radial coordinate 
and momentum corresponding to the classical turning 
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point respectively. Transforming the variable of integra- 
tion from p to e (using the relation /2m — t~ ^(t")) 
we get 

G(e) = y r^dr y v/2m(e - V^(r))de. (25) 

Thus the density of states is 



Jo 



rV2™ (e - Ar4)dr. (26) 



Substituting = x we can transform the integral into a 
form which can be evaluated analytically [l^ 



5(e) = 



/l3 



4V27r2m3/2r(3/2)r(3/4)e5/4 

;i3A3/4r(9/4) " 



(28) 



Using this expression for the density of states in eq. ^ 
we get 



iV4/9 



[r(9/4)C(9/4)]^/^ 



/i3A3/4r(9/4) 



r(3/2)r(3/4)4V27r2m3/2 



4/9 



(29) 

Comparing with cq. ([8]) we can obtain the ratio of Tc in 
the two cases, with and without the cross terms, for the 
3D isotropic quartic potential. It is found that Tc with 
the cross terms is 1.2 times higher. This rise in Tc can 
be attributed to the contribution from the cross terms 
which were neglected while deriving eq. ([5]). The reason 
for the difference is, when ^(e) = C^e""^ then Tc varies 
as Hence, the lower Tc in 3D quartic potentials 

without the cross terms is due to the higher value of Cq . 



B. 2D case 

In the 2D case, neglecting the cross terms, the potential 
is of the form 



V{x,y) = X{x^+y^). 



(30) 



Using the same approach as adopted in the 3D case, the 
total number of states available to the system is 



G{e) 
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^3/2 / ^i-'^'de^ / 6-^/^*2(31) 
16(ftw)3/^ Jq Jq 



3r(3/4)r(7/4) 



Thus the density of states 

0.9531ei/2 



?(e) 



(fia')3/2 



(32) 



(33) 



Substituting this expression of g(e) in eq. ([7]) we get 



kTc^ 



[C(3/2)r(3/2)0.9531 



i2/3 ■ 



(34) 



The corresponding expression when cross terms are con- 
sidered is fill 



kTr 



27r2TOr(3/2)C(3/2) 



2/3 



(35) 



Comparing eq. and eq. ([55]) we find that with the 
cross terms in eq. psp is approx. 1.12 times higher than 
Tc without the cross terms in eq. Thus the cross 

terms increase Tc in 2D as well as 3D case. 



CONCLUSIONS 



Our calculations show that Tc in the case of the 3D 
quartic potential trap is higher than that of the isotropic 
harmonic potential trap. This is due to the form of the 
density of states g{e), which varies as e^/^ and in 3D 
isotropic quartic and quadratic trapping potentials re- 
spectively. This implies lower density of states in quartic 
oscillator potential compared to isotropic harmonic oscil- 
lator potential. However, more interesting is the effect of 
the cross terms. In the 3D isotropic harmonic potential 
trap the cross terms are absent, which is not the case for 
the 3D quartic potential trap. The cross terms tend to 
decrease the density of states and raise Tc- These terms 
increase Tc by factor of 1.2 and 1.1 in the 3D and 2D quar- 
tic trap potentials respectively as compared to the case 
without the cross terms. Experimentally, in optical traps, 
the potentials without the cross terms are more appro- 
priate. We find that the effect of finite particle number 
is more pronounced in the 3D quartic potential when the 
number of particles is < 10^. The situation is reversed 
when the number of particles is > 10^. The cause of the 
reversal lies in the form of the fractional change ATc/Tc 
for the two potentials. The ratio of the fractional change 
between 3D isotropic quartic potential to harmonic po- 
tential is 3.51/iVi/9. It is « 1 for N « 10^ > 1 for 
N < lO'^ and < 1 for N > 10^ Thus, when N < 10^ the 
constant factor is dominant in eq. (j20p and is responsi- 
ble for the larger value of ATc/Tc in quartic potential. 
But when N > 10^ TV-^/g dominates and ATc/Tc of 
the quartic potential is lower than that of the harmonic 
potential. 
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